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Abstract. 

We report measurements of the evolution of lines, planes, and volumes in an 
intensely turbulent laboratory flow using high-speed particle tracking. We find that the 
classical characteristic time scale of an eddy at the initial scale of the object considered 
is the natural time scale for the subsequent evolution. The initial separation may only 
be neglected if this time scale is much smaller than the largest turbulence time scale, 
implying extremely high turbulence levels. 
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The transport of material by a carrier fluid is ubiquitous in both the natural world 
and in engineering applications [1]. When the carrier flow is turbulent, the dispersion 
of the transported substance can be very rapid. Turbulent flows are also extremely 
efficient at mixing [2], since their nonequilibrium nature drives the production of small 
scales and sharp gradients where diffusion can occur rapidly. To study both transport 
and mixing, it is natural to work in the Lagrangian framework where the fundamental 
objects are the trajectories of individual fluid elements [3]. 

The trajectory of a single Lagrangian particle is, in general, not sufficient for 
characterizing turbulent transport. Instead, knowledge of the collective motion of groups 
of particles is required [1,4]. The simplest such multiparticle quantity is the growth of 
the relative distance between a pair of particles [5]. We have previously measured 
this relative dispersion [6,7], finding excellent agreement with Batchelor's theoretical 
predictions [8]. Geometrically, the two particles in relative dispersion define a line, 
and therefore a particle pair only gives one-dimensional information about turbulent 
dispersion. An even deeper understanding of turbulent transport requires the study of 
higher-dimensional structures: groups of three particles, which define a plane, and of 
four particles, which define a volume. While such structures have been considered before 
in models and numerical simulations [9-1 1] and in low- Reynolds-number experiments 
[12,13], they have not been investigated at the high Reynolds numbers that are common 
in nature. 

Here, we present measurements of the shape dynamics of collections of Lagrangian 
particles in an intensely turbulent laboratory water flow. We consider lines (two 
particles), planes (three particles), and volumes (four particles). In all of these cases, 
the initial size of the object plays a strong role: the time scale determined by this initial 
size, which we denote by t , characterizes the experimentally observed evolution [14]. We 
find that two particles separate superdiffusively, as is well-known in turbulence, but that 
to divides two types of separation behavior. Triangles formed from three particles and 
volumes spanned by four particles also grow in time, but assume stationary shapes after 
to, with triangles evolving to a preferred set of internal angles and volumes flattening 
into nearly planar structures. 

Our measurements are made using optical particle tracking [15] in a swirling water 
flow between counter-rotating baffled disks, as described in detail elsewhere [7]. We 
characterize the strength of the turbulence with the Taylor-microscale Reynolds number 
R\ = y/lbu'L/ v, where u' is the root-mean-square velocity, L is the correlation length 
of the velocity field, and v is the kinematic viscosity; here, we report measurements for 
R\ as high as 815. Our polystyrene tracer particles are smaller than or comparable to 
the smallest scale of the turbulence, the Kolmogorov length scale rj = (z/ 3 /e) 1//4 , where e 
is the mean rate of energy dissipation per unit mass, for all Reynolds numbers studied, 
and faithfully follow the flow [16]. By modifying our tracking algorithms, we have been 
able to increase the lengths of measured trajectories significantly, allowing for the study 
of longer-time statistics [17]. 

Let us first consider one-dimensional shape changes by measuring the growth of 
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Figure 1. (a) The separation of two particles in time, compared with a modified 
Richardson-Obukhov law, for different initial separations, ranging from 1 mm to 5 
mm. The Reynolds number is fixed at R\ — 690, with a Kolmogorov scale of n = 30 
/im. We observe similar behavior at other Reynolds numbers, (b) The change of Cr 
with T^/to, shown for three Reynolds numbers. 



the separation R(t) between two particles. For this case, the well-known Richardson- 
Obukhov law [5, 18] predicts that 

(R 2 (t)) = get 3 (1) 

in the inertial range, i.e., rj <C R <C L and r v <^ t Tl, where Tl = (L 2 /e) 1//3 is 
the large-eddy turnover time. In Eq. (0Q) the dimensionless coefficient g, known as the 
Richardson constant, is expected to be universal and independent of initial separation 
i?o- It has been notoriously difficult, however, to observe conclusive evidence of this t 3 
scaling experimentally [5]. In our previous measurements of relative dispersion [6,7], we 
instead found that the initial separation of the pair R plays an important role, as first 
suggested by Batchelor [8]. He predicted that 

(SRiSRj = (ll/3)C 2 ( ej R ) 2/3 t 2 , (2) 

where 5H(t) = R(£) — R and C2 = 2.13 ± 0.22 is the scaling constant for the second- 
order Eulerian velocity structure function [19]. Batchelor additionally predicted that 
this scaling law should hold for t <C to, where 

to = (i? 2 A) 1/3 , (3) 

may be regarded as the lifetime of an eddy of scale R . 

It has been suggested that the failure to observe the Richardson-Obukhov law is 
due to the influence of particle pairs that separate anomalously slowly or quickly, so that 
they bring in non-inertial-range effects [20,21]. This will occur unless the inertial range is 
sufficiently wide and the effects of both the dissipation and integral scales are negligible, 
requiring very large Reynolds numbers. In addition, the finite measurement volume in 
experiments may introduce a bias against quickly separating particle pairs [22]. We 
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Figure 2. (a) Change of the energy of relative motion 5u 2 (t) — 5u 2 (Q) following 
the pairs, normalized by (ll/3)(ei?o) 2 ^ 3 , the theoretical average relative energy at the 
initial separation between the pairs, where 8u is the relative velocity between the pairs, 
(b) The rate of energy change dE/dt normalized by the turbulence dissipation rate, 
where E = (Su 2 )/2. Both (a) and (b) are from the R\ = 690 experiment, with a 
Kolmogorov scale of rj = 30 /im. The legends give the initial separations in mm. Wc 
observed similar behavior at other Reynolds numbers. 



have therefore also measured (R 2 ^ 3 (t)} — Rq , which is less affected by the finite- 
volume bias and may display scaling behavior at Reynolds numbers accessible in current 
experiments [6]. If the Richardson- Obukhov law holds, then 

((R 2/3 (t)} - RT)/RT = C R (t/t ), (t < t < T L , R < R < L) (4) 

where Cr should be a constant related to the Richardson constant g. The compensated 
plot ((.R 2//3 (£)) — Rq)/Rq /it /to) should collapse to a plateau in the inertial range, 
independent of initial separation Rq. As shown in Fig. H{a), plateaus, though short, 
do exist for t > f . We find, however, that the initial separation R again plays a 
role: to is the time scale of the transition to the (i? 2//3 ) ~ t scaling and the plateau 
values depend on Rq. These observations support our earlier argument that a very 
large separation between Tl and to (corresponding to a very large Reynolds number) 
is required to observe i?o-mdependent Richardson-Obukhov scaling [6], which is also 
supported by recent work using a stochastic model [22]. To quantify the effect, we plot 
in Fig. [U(b) the change of Cr with Tl /to at three different Reynolds numbers, where 
Cr is measured from the peak value in the plateau region of the compensated curves. 
The largest Ti/to in this plot are taken from Rq m 30r), the lower bound of the inertial 
range [23] . It can be seen that the measurements from three different Reynolds numbers 
collapse very well, and that the _R - m dependent region is not reached even at R\ = 815. 

The absence of Richardson universal scaling law even at R\ ~ 10 3 may be better 
understood from the evolution of 5u, the relative velocity between fluid particles pairs. 
Figure[2](a) shows the change of the energy of relative motion, 5u 2 {t) —5u 2 (0), following 
the pairs. Figure W(Jd) shows the rate of the change, dE/dt, following the pairs, where 
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Figure 3. The evolution of Lagrangian triangles at R\ = 815 (77 = 23 //m), 
characterized by (a) the mean-squared lengths of the triangle edges and (b) the mean 
triangle angles. Rq, the initial edge length, varies from 4 mm to 6 mm. As in the 
two-particle case, the data for different initial sizes collapse when time is scaled by to. 



E = Su 2 /2. Clearly, the energy of relative motion decreases initially before the final 
increase. As initial separation increases, the return to dE/dt > takes longer time. 
These observations are in good agreement with results Numerical simulations [24] , where 
the same change of the energy of relative motion was found for a "cloud" of tracer 
particles. This initial decrease of relative motion energy can be linked to the large 
correlation length of forcing in three-dimensional turbulence [4] . In light of Figure [21 
the relative motion between a pair of particles is first slowed down and it takes a time 
comparable to to for the relative motion to accelerate again. The further apart the 
particles the longer the relative period of slow down is. Only after the initial decreasing 
period, the Richardson regime can possibly be observed. 

While studying the two-particle case tells us how particles separate along lines in 
turbulence, groups of three particles can show richer dynamics since they form two- 
dimensional objects. We therefore measured the shape statistics of triangles formed 
from three particles. Here we report only the statistics of triangles that were initially 
nearly equilateral. We then measured the growth of the sides the triangles, analogous 
to the evolution of the mean-squared separation of two particles considered above, and 
the distortion of the triangles, information that is unavailable from pair separation 
measurements. 

Figure [3](a) shows the evolution of the mean-squared lengths of the triangles (Af ) 
(with Ai > A 2 > A 3 ) for three different initial triangle sizes. When time is scaled by t , 
defined now with the side length of the initially equilateral triangles, the data collapse 
for different initial sizes. Additionally, scaling in this fashion collapses data taken at 
different Reynolds numbers (not shown). We note that unlike the monotonic increase 
of (Af) and (A|), there is initially a slight decrease in (A3), an effect that is absent in 
the statistics of particle pairs. This result shows the importance of studying higher- 



Evolution of Geometric Structures in Intense Turbulence 



6 




Figure 4. The evolution of Lagrangian tetrads at R\ — 690 (77 = 30 /mi), 
characterized by (a) the mean eigenvalues of the inertia tensor scaled by the initial 
tetrad size and (b) the mean normalized eigenvalues for initial sizes ranging from 10 
mm to 20 mm. Once again, the data for different initial sizes collapse when time is 
scaled by to. 



dimensional structures in three-dimensional turbulence. We shall discuss this initial 
decrease in more detail when considering four-particle statistics below. 

In addition to a mean growth, three particles can assume nontrivial 2D 
configurations. To quantify these shapes, we show the evolution of the means of the three 
triangle angles (0$) in Fig. 13(b). Just as for the edge lengths, the angle data collapse 
when time is scaled by to- We also observe that our initially equilateral triangles become 
distorted over a time of order to to stationary, obtuse shapes, with angles of 0.567T, 0.277T, 
and 0.177T. These values appear to be independent of Reynolds number: indeed, our 
observations are in agreement with earlier results from 3D numerical simulations [10] and 
2D experiments [12], even though they were performed at significantly lower Reynolds 
number. 

The evolution of triangles has given us more insight into transport than the two- 
particle case, but for 3D turbulence, we must consider volumes for a full characterization. 
We therefore consider groups of four particles, the minimum number of points required 
to define a volume. Such groups form "tetrads," which have previously been used to 
construct a stochastic model of the coarse-grained velocity gradient tensor [9]. Following 
Chertkov et al. [9], we characterize the spatial arrangement of the particles with the 
eigenvalues g^ (with g\ > g 2 > #3) of the "inertia" tensor = X^X^j. The column 
vectors of the tensor Xij are defined based on the relative separation between the four 
points: 

Xi = (X 2 - Xl ) /y/2, 

X 2 = (2x 3 -x 2 -x 1 )/v / 6, (5) 
X 3 = (3x 4 - x 3 - x 2 - xi) / VV2, 
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(a) (b) 

Figure 5. The evolution of PDF of the shape factor following Lagrangian tetrads 
with initial size Rq = 20mm in turbulence. The legends show the time (in units of r^) 
at which the PDFs are measured. Similar changes are observed for tetrads with initial 
size of 10 and 15mm. (a) R x = 690 (77 = 30 ^m); (b) R x = 815 (77 = 23 /mi). 



where x ra (n = 1,2,3,4) is the position of the n th Lagrangian particle. The volume of 
the tetrad is given by V = {^/gig^gt) /3, and its radius of gyration is R 2 = g\ + g 2 + g% 
[10,11]. We also define normalized eigenvalues Jj = g%jR 2 g to describe the tetrad shape: 
I\ = I2 = I3 = 1/3 defines an isotropic tetrad, I3 = means that the four points are 
coplanar, and I 2 = I3 = means that they are collinear. As with the triangles studied 
above, we considered only tetrads that were initially nearly isotropic, with the length 
of all edges within 10% of a nominal size Rq. 

In Fig. H](a), we show the evolution of the tetrad eigenvalues scaled by the initial 
tetrad size. Data for different initial sizes collapse when time is scaled by to, as it would 
for different Reynolds numbers. Here, it is clear that while the two larger eigenvalues 
increase monotonically, the smallest eigenvalue decreases. This observation provides 
direct experimental evidence for earlier theoretical and numerical work [9] that in 3D 
turbulence the large-scale, coarse grained strain-rate tensor has, on average, two positive 
eigenvalues and one negative eigenvalue. The negative eigenvalue causes the initial 
compression in one direction, and consequently the decrease in (#3). The normalized 
eigenvalues, shown in Fig. H^b), also collapse for different sizes when scaled by t . We 
also find that the tetrads develop towards a nearly planar configuration after a time of 
order to- Within the experimental resolution, we determine the stationary values of the 
normalized eigenvalues to be (I 2 ) = 0.25 and (J 3 ) = 0.06 (note that (Ii) = 1 — (h) — (I3)), 
similar to our earlier results where we did not select only the initially nearly isotropic 
tetrads [25]. We note that we did not observe a Richardson- like ^ ~ t 3 for any Reynolds 
number or initial size accessible. 

Our results are qualitatively similar to the recent numerical simulations of 
Lagrangian tetrads in Ref. [11], performed at R\ = 280, where stationary values of 
(I2) = 0.16 and (I3) = 0.02 were reported, indicating an even stronger tendency towards 
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Figure 6. The evolution of the averaged tetrad volume (V), normalized by Rq, 
following Lagrangian tetrads with inertial range sizes in turbulence, (a) R\ = 690 
(jj = 30 fmx)\ (b) R x = 815 (77 = 23 fmi). 



planar tetrads than in our data. In Fig. [5j we show the evolution of the probability 
density function (PDF) of the smallest shape factor (J 3 ) following Lagrangian tetrads. 
Again, the evolution of the PDFs are consistent with simulation results in Ref. [11] 
for tetrads with initial separation in the dissipative range. It can be seen that after 
60 ~ 70t v (or approximately t ), the PDFs are still slowly evolving. Beyond that time, 
the tetrads, on average, have been swept away from their initial positions by half of 
the size of the measurement volume and disappear from our observation. Therefore, our 
measurement of (ij) might be affected by the finite measurement volume. In addition, it 
is more probable for a vertex of an extremely deformed tetrad to move out of the finite 
measurement volume. There is hence a potential bias towards smaller deformation in 
our experiments. 

A possible reason for the slight difference between the simulation and experimental 
results could be the higher potentiality for a vertex of an extremely deformed tetrad 
to move out of the finite measurement volume of our experiments, at which point the 
tetrad would no longer be detectable. There is therefore a potential bias towards smaller 
deformation in our experiments. 

The initial size of the tetrads studied in Ref. [11] are in the dissipation range 
(Rq ~ 77), whereas our smallest initial separation was Ro ~ 300?7, well into the inertial 
range. In Ref. [11], the initial decrease in #3, similar to that shown in Fig. Ufa), was 
attributed to dissipative effects, which were also presumed to be responsible for the lack 
of observed Richardson scaling. Since we observed very similar behavior with Rq well 
in the inertial range, however, we attribute the lack of t 3 scaling to a small separation 
between to and Tl, just as for the case of pair separation. The correspondence between 
our inertial-range results and the smaller-scale numerical results suggests additionally 
that the statistics of the coarse-grained velocity gradient tensor [9] may be very similar to 
those of the true velocity gradient. Further study of these similarities will be an excellent 
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opportunity to bring together complementary experiments and numerical simulations. 

Figure [6] shows the change of averaged volume following the tetrads with different 
sizes. The tetrads volume remains constant initially for time t < O.lio- Then decreases 
in volume before the final increase. This decrease was not observed in previous numerical 
work [10] and experimental measurements [13]. A possible reason is that the tetrad sizes 
are large enough to feel the effect of the mean strain. 

In summary, we have measured multiparticle Lagrangian statistics in an intensely 
turbulent laboratory flow. We studied the evolution of lines, planes, and volumes, 
each parameterized by the minimum number of Lagrangian points. We find that, in 
three-dimensional turbulence, material volumes tend to flatten into planar shapes, in 
agreement with previous numerical and experimental studies at much lower Reynolds 
numbers. Our results clearly indicate that the initial particle separation is an important 
parameter in each of these cases and must be included in models. Only when the ratio 
of the largest turbulence time scale and the time scale based on the initial size of the 
object is very large, implying a very large Reynolds number, may the initial size be 
safely neglected. 
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